Transient nanobubbles in short-time electrolysis 
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Water electrolysis in a microsystem is observed and analyzed on a short-time scale ~ 10 fis. Very 
unusual properties of the process are stressed. An extremely high current density is observed because 
the process is not limited by the diffusion of electroactive species. The high current is accompanied 
by a high relative supersaturation S > 1000 that results in homogeneous nucleation of bubbles. 
On the short-time scale only nanobubbles can be formed. These nanobubbles densely cover the 
electrodes and aggregate at a later time to microbubbles. The effect is significantly intensified with 
a small increase of temperature. Application of alternating polarity voltage pulses produces bubbles 
containing a mixture of hydrogen and oxygen. Spontaneous reaction between gases is observed 
for stoichiometric bubbles with the size smallaer than 150 nm. Such bubbles disintegrate violently 
affecting the surface of electrodes. 

PACS numbers: 47.55.D-, 64.60.qj, 82.33.Vx 
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I. INTRODUCTION 

More than 200 years gone since people observed wa- 
ter electrolysis for the first time (see [l|, Q and references 
therein). Since then the process is widely used in dif- 
ferent applications including hydrogen production 0, 0] ■ 
Nowadays electrolysis of water find also applications in 
different kind of microdevices such as actuators [HI], 
pumps [1, [l(| , and others [HI ■ Fast performance of 
microsystems requires processing on a short-time scale. 
Surprisingly enough, most of the researches on water elec- 
trolysis dealt with macro or microscopic systems on a 
time scale longer than milliseconds [HEl. A few papers, 
where shorter (microseconds) voltage pulses were con- 
sidered, discussed increase of the overall efficiency 1141 
or optimization of the performance of electrolysis [15j . 
but did not observe the process on the microsecond time 
scale. 

In normal long-time electrolysis (long-time means here 
the time scale r > 1 ms) gas bubbles with the diam- 
eters in the range 10 — 1000 jum are observed [HI, fl6| . 
The maximal current density that can be reached in this 
case is around 1 A-cm ~ 2 Il7l| and the maximal relative 
supersaturation ~ 100 [l8|]. These values cannot be in- 
creased further because bubbles cover all the surface of 
the electrode making it inactive [13, OH, Hq| • 

In the short-time electrolysis, say t = 1 — 100 /is, one 
could expect unusual properties of the process. It be- 
comes obvious if we consider gas diffusion into the liq- 
uid above the electrode. For the time r a significant 
supersaturation is possible at the distance from the elec- 
trode I ~ \J Dt, where D is the diffusion coefficient of 
the formed gas in the surrounding liquid. Taking as typ- 
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ical values D ~ 10~ 9 m 2 /s and t ~ 10 fis one finds 
I ~ 100 nm. It means that on this timescale only 
nanobubbles can be formed. This is indeed the case as we 
demonstrated recently [2l|. Nanobubbles of 200 — 300 nm 
in diameter were observed stroboscopically above mi- 
croelectrodes. These nanobubbles were formed for the 
first 20 — 50 ^s then aggregated to microbubbles after 
100 — 300 /its, but the following dynamics was very slow. 

A number of observations indicate that smaller 
nanobubbles can exist at shorter times r < 10 /us. For 
example, it was demonstrated that the reaction between 
hydrogen and oxygen can be ignited spontaneously in 
nanobubbles with the size smaller than 150 nm [2l|. Dif- 
ferent manifestations of this reaction show that the reac- 
tion persists for nanobubbles as small as 50 nm or even 
smaller. 

In this paper we describe formation of nanobubble on 
the short-time scale r = 1 — 100 /is. It is demonstrated 
that very high current densities and supersaturations 
can exist accompanied by the homogeneous nucleation of 
nanobubbles. The bubbles containing hydrogen or oxy- 
gen grow and aggregate forming micro-sized bubbles on 
the timescale r > 100 /is. The nanobubbles contain- 
ing mixture of the gases disappear as the result of the 
reaction between gases. These are transient nanobub- 
bles, which exist only a short time and difficult for direct 
observation. In contrast with the first report [2l| we pro- 
vide here detailed description of the transient nanobub- 
bles and pay less attention to the chemical reaction inside 
of these bubbles. 

The nanobubbles discussed in this paper differ from 
the surface nanobubble, which attracted recently a con- 
siderable attention. The long-lived surface nanobubbles 
were observed with different methods [22T4241 ] . Observa- 
tion with atomic force microscopes (AFM) became es- 
pecially popular (see [25[ and references therein) due to 
its simple realization. These nanobubbles are widely dis- 
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cussed in the literature because of their mysterious sta- 
bility. High Laplace pressure drives diffusion of gas out 
of the bubble resulting in a dissolution time of ~ 10 fxs 
but not days as observed (see [26| for a review). Precise 
reason for this stability is not clear yet but in a promising 
model [27j of " dynamical equilibrium" the gas flux going 
out of the bubble through the spherical cap is balanced 
with a gaseous influx at the contact line. 

A number of methods to generate the surface nanobub- 
bles are in use including solvent exchange, tempera- 
ture change for liquid or substrate, and liquid pressur- 
izing [Hj]. The electrochemical formation of the surface 
nanobubbles were observed in water electrolysis for hy- 
drogen Q and for both hydrogen and oxygen [30| on 
highly orientated pyrolytic graphite. The electrochemi- 
cal process is considered as a flexible way to control the 
bubble production (density and size) [311 ] . In contrast 
with the electrochemical process considered in this paper 
the electrochemical production of the surface nanobub- 
bles had much smaller current density (~ 10~ 4 A/cm 2 ) 
and much longer time scale (~ 10 s). 

This paper is organized as follows. In Sec. [TT] we de- 
scribe the system for observation of the short-time elec- 
trolysis, application of the voltage pulses, and the cur- 
rent response. In Sec. Mil the main results obtained with 
the stroboscopic observations are described. Electrolysis 
with alternating sign pulses is described in Sec. II VI where 
we present the results obtained with a stroboscop and vi- 
brometer. We explain also modification of the electrode 
surface due to the process. In Sec. [V] the evolution of 
the transient nanobubbles is summarized and we discuss 
their possible applications. Our conclusions are collected 
in the last section. 



II. SHORT-TIME ELECTROLYSIS 

To analyze electrolysis on a short-time scale in a mi- 
croscopic system special microchips were prepared. One 
of the chips is shown in Fig. [T] For fabrication we used 
wafers made of borofloat glass BF-33. First, Ti sublayer 
was deposited on glass (10 nm) for better adhesion, then 
the 100 nm thick metallic layer was sputtered. We used 
the wafers with different metals on top such as Pt, Pd, 
Au, or W. The metals were patterned and covered litho- 
graphically with the insulating resist SU8 where it was 
necessary. Different thicknesses of SU8 layer of 3 and 
90 /xm were in use. The fabricated wafers were diced on 
separate chips of size 12 x 12 mm 2 . Each chip contained 
16 pairs of the microelectrodes of different shapes and 
sizes and the contact pads to address a specific pair elec- 
trically. The patterned SU8 layer formed walls of open 
channels for a liquid. A chip filled in with the electrolyte 
was covered by a thing (30 /im) glass plate for observa- 
tion of the electrolysis. 

As the solution for electrolysis we used 15 g of Na 2 S0 4 
dissolved in 100 ml of deionized water. Also we have 
used the solutions of NaCl and KI in water in similar 




FIG. 1: (a) Microchip used for investigation of electrolysis 
on short-time scales, (b) Pair of microelectrodes that can be 
addressed separately. 

concentrations. 

The process was observed with a homemade strobo- 
scopic system [32| . In this system the light source and the 
electrical pulses applied to the electrodes were controlled 
with the time resolution 0.1 /is. As a light source two 
powerful green LEDs LUXEON® III were used. With 
this system we have been able to observe the electrodes in 
an optical microscope with the illumination time longer 
than 5 /js. 



A. Voltage 

For a chosen pair of the electrodes one was always 
grounded and the other one (working electrode) could be 
at a negative or positive potential. Hydrogen or oxygen 
bubbles were produced above the working electrode. We 
could apply voltage as a pulse of a given length from 1 
to 10 5 /is. The pulse could be repeated with a frequency 
/. At a prescribed moment of time a short flash of light 
> 5 /is was on and an optical microscope image of the 
electrodes was made during the flash. The applied volt- 
age and the resulting current through the electrodes were 
recorded. Alternatively keeping one electrode grounded 
we applied voltage pulses of alternating polarity to the 
working electrode repeated with some frequency. In this 
case bubbles containing both gases H 2 and 2 could be 
formed [2lj |. 

The voltage needed to start the electrolysis was al- 
ways larger than that for macrosystems, Uq = 2.5 — 4 V. 
The high threshold voltage was already observed in mi- 
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crosystems @, H, [HI , but the reason was not explained. 
It has to be stressed that in the microsystems a planar 
geometry is used when both electrodes are in the same 
plane. In this case the ohmic losses (resistance between 
the electrodes) are considerable even for well conducting 
electrolytes. This is in contrast with the macrosystems, 
where these losses are usually negligible. For one pair of 
electrodes we measured the resistance to be R = 2.3 kSl 
but it varies to some degree from pair to pair due to 
different geometry. The working current in our system 
I ~ 1 mA was quite large. Then the ohmic polarization 
is estimated as a few volts (2 — 3 V). 

Conventionally we defined the threshold for the elec- 
trolysis applying different voltages for 5 ms and observing 
the gas formation. At U = Uo the first gas was visible 
in the system. For example, for the same pair of elec- 
trodes, for which the resistance was measured, we found 
Uo = 3.2 V. In some cases Uq was as large as 6 V or so, 
but we have found out that this was due to the potential 
drop on the electrode surface. Treatment of such chips 
in 02-plasma significantly reduced Uq. Probably a film 
(hydrocarbons) was formed on top of the electrodes dur- 
ing the fabrication process. It should be stressed that the 
absolute values of the potential has to be taken with care 
because we deal with the two-electrode system and there 
is uncertainty in the potential of the working electrode 
in contrast with the three-electrode system. 



B. Current 

The current response to square voltage pulses applied 
to the working electrode is shown Fig. [21 The response 
is in accordance with the expectations. One can sepa- 
rate two components in the current. The first one is the 
faradaic current Ip that does not depend on time t for 
square voltage pulses. The second component is respon- 
sible for the charging-discharging of the double layer on 
the electrode surface and is not related to the electro- 
chemical reaction. It behaves as Iie _t / Tc , where I\ is a 
constant and t c is the relaxation time related to the ca- 
pacitance of the double layer. As the result we fit the 
current response on the square voltage pulse by a func- 
tion of three parameters: 

I{t)=I F + I 1 e- t /^. (1) 

It has to be stressed that typically in macro and mi- 
crosystems considered on the long-time scale the faradaic 
current Ip depends on the time because it is diffusion 
limited. For parallel electrodes the current is propor- 
tional to \jTo/t in accordance with the Cottrell equation 
[33| . where T£> is the diffusion time for the electroactive 
species in the solution. In our case the diffusion time is 
too long to play any role. Roughly it can be estimated 
as t d ~ L 2 /D a ~ 100 ms, where L ~ 10 fim is a typical 
size of our electrodes and D a ~ 10~ 9 m 2 /s is a typical 
diffusion coefficient of the electroactive species in water. 
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FIG. 2: (top) Current response to the positive square voltage 
pulses of the amplitude U = 5.2 V and frequency / = 40 kHz. 
(bottom) Response to the voltage pulses of alternating po- 
larity U = 6 V, / = 40 kHz. The insets show the currents 
from the main graph for half of one period and the fits of the 
currents with Eq. \T\ 

We used Eq. ([T]) to fit the current response as shown 
in Fig. [2J The quality of the fit can be seen in the insets. 
The faradaic current can be well approximated by a con- 
stant in the investigated time domain t < 100 /is. For 
alternating pulses the positive I p (t) and negative I m {t) 
halves of the pulse are symmetric as shown in the bot- 
tom inset. For single polarity pulses the relaxation time 
r c is larger for the negative half of the pulse (4.9 against 
1.3 (is) reflecting dependence of the double layer capac- 
itance on the potential [33| . For alternating polarity 
pulses in Fig. [5] we have found from the fit r c = 4.0 fis. 
With the frequency increase the relaxation time decreases 
slowly indicating that not all adsorption states can be 
filled for a short time. It has to be noted that the Fara- 
day component of the current persists up to the shortest 
pulses of 1 fj,s (/ = 500 kHz), which we investigated. 
Presence of this component signals formation of gases on 
this shortest time scale that was directly observed with 
the stroboscope. 

The faradaic current reaches rather large values. For 
the case shown in Fig. [2]it is Ip = 0.80 ±0.07 mA for the 
single polarity pulses and Ip — 1.08 ± 0.03 mA for the 
alternating polarity voltage pulses. Keeping in mind that 
the open area of the electrode was around 1260 ^m 2 one 
finds that the current density was jp ~ 80 A/cm 2 . This 
is the average value of the current density but locally it 
can be 5 times larger (see Sec. IIII Ap . This value is much 
larger than the maximal value ~ 1 A/cm 2 [ijj that can 
be reached in the long-time electrolysis. 

The most important conclusions that can be drawn 
from the current-voltage behavior of the microsystem on 
the short-time scale are the following, (i) The electrolysis 
persists normally up to the shortest investigated time 
1 /is. (ii) The current density is at least two orders of 
magnitude larger than was previously observed in the 



macrosystems. 



III. STROBOSCOPIC OBSERVATIONS 

The stroboscopic technique gives images that were 
made for different runs. To make conclusions from these 
images the system has to behave reproducibly from run 
to run. For our system this condition was fulfilled. The 
images made at the same moment for different runs were 
hardly distinguishable on the short-time scale and varied 
slightly on the long-time scale > 1 ms. In the latter case 
the positions and sizes of separate bubbles were slightly 
different but qualitatively the pictures were the same. 
For good visibility we used typically a 10 fxs long flash 
of green light that was delayed for different times from 
the beginning of the voltage pulse. The camera collected 
the light from the microscope during the flash. The es- 
timated time uncertainty in 5 /is is reasonably small al- 
though comparable with the time scale. The gas pro- 
duced by series of pulses in most cases disappeared dif- 
fusively in less than 1 s. Sometimes when the amount of 
produced gas was large a stable (pinned) bubble formed. 
In this case we have to refill the chip for further use. Due 
to evaporation we have been able to operate with a filled 
chip during 15 min. 

Figure [3] shows the electrodes in different moments of 
time. A positive or negative voltage pulse of 50 /zs long 
was applied to the electrodes. The left image shows both 
electrodes 20 /js after beginning of the positive voltage 
pulse. Oxygen and hydrogen are formed above the work- 
ing electrode and grounded electrode, respectively. All 
the other images show the working electrode with hy- 
drogen in different moments during the negative voltage 
pulse. 

Indeed, the periphery of the electrode is a preferable 
place for bubble nucleation, but one can see some haze in 
the inner part of the electrode. With our microscope we 
can reliably resolve bubbles with a diameter of 500 nm. 
The bubbles forming the haze are smaller. We are able 
to see them because they strongly scatter the light. The 
diameter of these scattering centers is estimated as d > 
X/tt s» 170 nm [34| where A = 520 nm is the wavelength 
of the used LED. These scattering centers fill densely the 
inner part of the electrode and grow with time. 



A. Supersaturation 

Let us make simple estimates. For the time r = 20 fis 
the gas molecules can diffuse in the vertical direction off 
the electrode on the distance 



Ihz ~ V d h 2 t ~ 300 nm, lo 2 « 200 nm, (2) 

where D H2 =4.5xl0~ 9 m 2 /sandL»o 2 =2.0xl0" 9 m 2 /s 
are the diffusion coefficients of hydrogen and oxygen in 
water. These values of the diffusion length have to be 
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FIG. 3: (left panel) Positive potential is applied to the work- 
ing electrode (bottom). The photo was made 20 fis after 
switching on the voltage. Six images on the right show the 
working electrode in different moments of time when negative 
voltage is applied to it. The indicated time is taken in the 
middle of the flash. 

considered as the upper limits on the bubbles that can 
be formed during the time r. Therefore, it is natural 
that we can see only nanobubbles during the first 20 fis 
of electrolysis. 

The relative supersaturation for H2 or O2 molecules 
can be found as follows. The concentration of i-th gas in 
the diffusion layer is rii = Afi/hA, where A is the area of 
the working electrode and Afi is the number of molecules 
produced by the electrolysis. This number is related to 
the Faraday current as Nh 2 = IpT/2\e\ (for H 2 ), where e 
is the electron charge. The relative supersaturation is the 
ratio of the gas concentration to the saturated concentra- 
tion nf \ which at normal conditions is 4.7 x 10 17 cm~ 3 
or 7.7 x 10 17 cm~ 3 for hydrogen or oxygen, respectively. 
Then we find for the supersaturation 



2 \ e \ n H 2 ^\ e \ n 2 

where jp — Ip /A is the current density. 

The faradaic current recorded together with the im- 
ages in Fig. [3] was If ~ 1-5 roA. However, it is not 
straightforward to estimate the current density because 
it is highly nonhomogeneous for the planar electrodes. 
We already analyzed the current distribution over the 
electrode (see [2l| Supplemental) extracting the informa- 
tion from the observed "wear" of the electrodes. It was 
found that the current density over the electrode surface 
can be presented as jp(x,y) = j m f(x,y), where x and 
y are the coordinates of a point on the electrode and j m 
is the maximal current density. The function f(x, y) has 
its average value over the electrode f(x,y) « 0.2. The 
area of the electrode was A ss 1260 /zm 2 and we find 
j m rj 600 A/cm 2 . This is an extremely high value in 
comparison with usual electrolysis and it is realized in 
the area of active bubble formation as can be seen in 
Fig. |31 The supersaturation in this area for the time 
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r = 20 /is is estimated as 

S% 2 ps 2700, Sg; ps 1200. (4) 
These are incredibly large values. 

B. Homogeneous bubbling 

Observation of the haze far from the electrode periph- 
ery is a signal of homogeneous nucleation of small bub- 
bles. Our estimates of the supersaturation (Q} also sup- 
port this idea. In the classical theory of homogeneous 
nucleation [H, [36| the probability to create a bubble of 
the critical size r c is given by the exponent 

f 47rrg 7 ) 2 7 

w ~ exp r^r)' rc = AP' (5) 

where 7 is the surface tension of the liquid and AP is the 
pressure difference between gas and liquid. The pressure 
in gas P g must support the high supersaturation in liquid. 
It is defined as P g = SP S , where P s is the saturated 
pressure for a gas at normal conditions. For both H2 and 
O2 gases P a ps 1 atm and we find APi ps P a (Si — 1) for 
the i-th gas. 

For small supersaturation S — 1 ~ 1 the critical size 
r c ~ 1 /zm is large and the exponent in ([5]) is incredibly 
small. For high supersaturation S ~ 1000 the critical 
bubble r c ~ 1 nm is small and the probability to cre- 
ate a critical bubble, which is able to grow, becomes ap- 
preciable. It is hardly reasonable to make more detailed 
estimates. This is because a number of important correc- 
tions have to be included in the classical nucleation the- 
ory [35ll37ll38| . It has to be noted that according to ([3]) Si 
increases with time as \fr. Therefore, starting from some 
moment after beginning of the electrolysis homogeneous 
bubble formation is expected. With the temperature in- 
crease as one can see from ([5]) the homogeneous bubbling 
must start earlier and due to the exponential dependence 
this tendency must be well pronounced. This prediction 
is very easy to check experimentally doing water electrol- 
ysis at elevated temperatures. 

For this experiment we heated the chip by a flat resis- 
tor and controlled the temperature by a thermocouple. 
A 20 (is long negative voltage pulse was applied to the 
working electrode. Images shown in Fig. @]were taken at 
the moments t = 17 fis (top row) just before the pulse end 
and at t = 27 /is (bottom row) 7 /is after the end of the 
pulse. Some chemical remnants are clearly visible on the 
electrode surface. They appeared due to imperfections 
in the fabrication process and visualize the area with the 
reduced ability for gas production. The temperature in- 
crease resulted in the increase of the Faraday current 
from I F ps 1.3 mA at T = 23° C up to I F ps 2.5 mA at 
T = 53° C. This increase is expected due to the increase 
of the ions diffusion coefficients. However, the increase 
of the bubble visibility and speed up of the dynamics are 
much more prominent. 




FIG. 4: Images of the electrode at t — 17 [is (top row) and 
at t = 27 /is (bottom row) taken at different temperatures. 
A negative voltage pulse of 20 /is long was the applied to the 
electrode. Hydrogen nano and microbubbles are visible. 

Increase of the solution conductivity (or equivalently 
current) with temperature was checked independently 
by measuring the conductivity of the bulk solution with 
Mcttlcr Toledo ScvcnMulti conductivity meter. We have 
found that the conductivity can be well described by the 
linear behavior a(T) = ao [1 + a(T — To)], where do is 
the conductivity at temperature To. The thermal coeffi- 
cient was found to be a — 0.024 K _1 . This value is in a 
good agreement with the value 0.028 K _1 found from the 
current increase on the chip and with the values reported 
in different literature sources. 

Comparing the images in the top row of Fig. 2] one can 
see that much more nanobubbles becomes visible with the 
temperature increase. On the other hand, the gas volume 
at T = 53° C must be approximately twice larger than 
at T = 23° C as follows from the current increase. It 
has the following explanation. At higher temperature 
the nucleation of critical bubbles happened earlier and 
to the moment t = 17 /is much more bubbles had enough 
time to grow to a visible size. At T = 23° C most of the 
bubbles arc small enough to be invisible. 

It becomes even more clear when one compares the 
top and bottom rows. For example, the bottom image 
at 23° C contains much more visible gas than the top 
image. However, when the pulse is over there is no gas 
production anymore and one would expect the equal or 
the opposite relation between the gas volumes. There 
is a very simple and natural explanation for this phe- 
nomenon. The voltage pulse produces a certain amount 
of gas. This gas can exist in three forms: molecules dis- 
solved in the liquid, nanobubbles larger than the critical 
size but smaller than needed to be visible with the strobo- 
scope, and visible nano or microbubbles. If homogeneous 
nucleation of bubbles happens, then at a certain moment 
a significant volume of gas is collected in the invisible 
nanobubbles. These nanobubbles grow with time and 
become visible with a delay producing the phenomenon 



6 





m~ 

&* 2* 

?> f> 




FIG. 5: Long-time dynamics of bubbles. Hydrogen (top row) 
or oxygen (bottom row) bubbles were produced by a voltage 
pulse of 100 (is long. The images in each row were taken at 
different moments of time. 

of gas appearing from nothing. 

We consider Fig. |4]and explanations above as a strong 
support of the idea that in the short-time electrolysis we 
observe homogeneous formation of bubbles. 



IV. ALTERNATING SIGN PULSES 

Up to now we discussed only hydrogen and oxygen 
bubbles that are formed on different electrodes when neg- 
ative or positive voltage is applied to the working elec- 
trode. However, there is a simple but effective procedure 
that allows formation of bubbles containing mixture of 
H2 and O2. Moreover, the gas composition can be con- 
trolled electrically. The idea is to use pulses of alternating 
polarity 2lJ. In this case hydrogen and oxygen will be 
produced locally above the same spot on the electrode. 
If the pulses are short enough to neglect the gas redistri- 
bution due to diffusion, then one can expect formation 
of bubbles containing both gases. 

For homogeneous nucleation formation of a bubble 
containing both gases is favorable in comparison with 
two bubbles containing different gases. This is because 
even in the regime of homogeneous nucleation there is an 
energetic barrier to form the critical nuclei. For the sec- 
ond gas it is easier to diffuse inside of the existing bubble 
than to overcome the barrier and form a new nuclei. 



Reaction inside of nanobubbles 



C. Transition to long-time dynamics 

Let us consider now how transition to the usual long- 
time electrolysis happens. A negative or positive voltage 
pulse of 100 (is long was applied to the working electrode. 
The electrode was observed in different moments of time 
when the pulse was already switched off. Some results 
are shown in Fig. [5j 

It was already clear from Fig. 2] that nanobubbles 
growing from nucleus start to coalesce in some moment 
of time and form microbubbles (see bottom row in Fig. 
2|). Fast aggregation becomes possible due to high den- 
sity of nanobubbles. Now from Fig. [5] one can see that 
significant changes happen between the images taken at 
100 and 300 (is. The following dynamics becomes very 
slow. This is especially clear for the images at 0.5 ms 
and 5 ms. It is expected that the long-time dynamics is 
controlled by the slow diffusion and slow coalescence. 



The main conclusion that one can make from the 
short-time stroboscopic observations is that the nanobub- 
bles are densely produced in the homogeneous nucleation 
regime. They exist for some time in the "invisible" form 
when they are small enough to scatter the light and are 
invisible for the stroboscope. At later moments they be- 
come visible but exist as separate nanobubbles. Later 
they start to aggregate to form microbubbles that are 
usually observed in the electrolysis. Extremely high cur- 
rent densities and relative supersaturations are observed. 



While the frequency of alternating pulses is not very 
high the process proceeds similar to that for single po- 
larity pulses. However, when frequency is higher than 
/ = 20 kHz (at room temperature) the visible gas sud- 
denly disappears. This can be clearly seen in Fig. [5] 
The figure shows gas developed in the system after 1 ms 
of electrolysis. Negative potential pulses applied to the 
working electrode (bottom in each image) result in H2 
formation above this electrode at different frequencies 
(right column). No significant dependence on the fre- 
quency is observed. Oxygen formed above the working 
electrode for positive pulses also does not show frequency 
dependence (middle column) . The left column shows the 
situation for the alternating pulses applied to the work- 
ing electrode. In this case we observe strong dependence 
on frequency. The alternating pulses repeated with the 
frequency / = 20 kHz produce significant amount of gas. 
One can see that this gas is collected in bubbles which 
differ from that in the middle and right columns. Already 
at / = 50 kHz only small volume of gas is observed. At 
/ = 100 kHz practically no gas is visible. 

Where the gas has disappeared? The Faraday com- 
ponent of the current depends on frequency very weakly 
if depends at all. Because the pulses of single polarity 
produce the gas at all frequencies it is difficult to assume 
that the Faraday's law breaks at higher frequencies. Ex- 
perimentally we established that the gas disappearance 
is related to the stoichiometric production of hydrogen 
and oxygen. There is no visible gas production while 
positive and negative halves of the pulse have the same 
duration. In this case two H2 molecules are produced for 
one O2 molecule. If we change the relative duration of 
the positive and negative parts of the pulse then the sto- 
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FIG. 6: Gas production with the negative pulses (right col- 
umn), the positive pulses (middle column), and the alternat- 
ing polarity pulses (left column) at different frequencies. For 
the alternating pulses the visible gas production disappears 
suddenly at higher frequencies. 

ichiometric balance will be broken. In this case the gas 
reappeared in the system as shown in Fig. [7] In this fig- 
ure D is the duty cycle of the pulse that is the fraction of 
time when voltage on the working electrode is negative. 

It seems reasonable to assume that the gas disappears 
due to a chemical reaction between gases happening in 
the solution. The reaction is hardly possible between sep- 
arate dissolved gas molecules. For high supersaturation 
the molecules exist in the solution only a short time be- 
fore entering nanobubbles as one can see in Fig. [3] and 
21 It is possible that the reaction proceeds in nanobub- 
bles containing stoichiometric mixture of H2 and 02- Of 
course, the bulk mixture of gases at room temperature 
and normal pressure will not react, but the gases con- 
fined in a small volume could behave differently. It is 
known that below a certain size, material properties can 
change drastically [39(. For example, the surface tension 




FIG. 7: Visible gas in the system after 200 /is of processing at 
/ = 100 kHz for different duty cycles D of alternating polarity 
pulses. When the working electrode is equal time at positive 
and negative potential (D — 0.5) no gas is visible. If each 
pulse is more negative (D > 0.5), then excessive hydrogen is 
formed. For D < 0.5 excessive oxygen appears. In all cases 
the deviation of D from 0.5 results in reappearance of the 
visible gas. 

can support metastable phases of nanocrystals [!(| that 
exist only at high pressure for bulk materials. In liq- 
uids the surface tension results in significant pressure in- 
side of nanobubbles. This pressure can shift the chemical 
equilibrium [4l| and it is known also that fast dynami- 
cal processe s p l ay role for ignition of the reaction on the 
macroscale p3 - |44| . Nevertheless, the precise mechanism 
of the reaction is not clear at this moment. 

We investigated also the effect of gas disappearance 
in different solutions such as Na2SC>4, NaCl, and KI 
dissolved in deionized water to similar concentrations 
(w 1 M). The results of alternating pulse electrolysis in 
these solutions are shown Fig. [5] In the case of NaCl 
hydrogen and chlorine are produced. These gases can re- 
act with each other in the exothermic reaction with the 
enthalpy AH = — 92 kJ/mol considerably smaller than 
that for hydrogen and oxygen AH = —242 kJ/mol. One 
can see that the amount of visible gas decreases with the 
frequency increase but does not disappears completely. It 
means that the reaction starts spontaneously in smaller 
bubbles than for Na2S04 solution. In the case of KI solu- 
tion hydrogen and iodine are produced. The reaction be- 
tween these molecules is endothermic, AH — 25 kJ/mol, 
and cannot proceed spontaneously. Only weak depen- 
dence on / is observed that is rather related to a small 
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FIG. 8: Alternating pulse electrolysis in different solutions: 
Na2SC>4, NaCl, and KI dissolved in water. Each column cor- 
responds to one of the solutions labeled on the top. Gases 
formed in the electrolysis are indicated in the bottom of the 
columns. Each row corresponds to the same frequency. The 
effect of gas disappearance clearly visible for Na2S04 solution 
becomes less pronounced for NaCl and disappears at all for 
KI solution. 

variation of the current. 

Below we provide some additional experimental infor- 
mation that supports the hypothesis of chemical reaction 
inside of nanobubbles. However, the reaction is not the 
main topic of this paper; more information can be found 
in El. 



B. Observations with vibrometer 

Observation of gas with a vibrometer is complimentary 
to observation with the stroboscope. The stroboscopic 
picture becomes visible when the bubbles significantly 



scatter the light. The nanobubbles with the size d < X/n 
are invisible for the stroboscope. On the other hand, 
strong scattering of the laser beam makes the interfer- 
ometric (vibrometer) measurements impossible. There- 
fore, the vibrometer is effective for very small bubbles 
d < X/ir where the stroboscope fails. 

The vibrometer is sensitive to variations of the optical 
path Ad(t) of the laser beam with time. The optical path 
is defined by the change of the refractive index An due 
to the presence of gas in the solution. If we assume that 
within the laser beam the lateral distribution of gas is 
homogeneous then 



Ad = J dz [Am(z) + An 2 (z)} 



(6) 



where z is the vertical coordinate counting from the elec- 
trode surface, and Arii(z) is the change of the refraction 
index due to gas i. Note that Ad does not change if the 
gas molecules redistribute along z but their number stays 
the same. 

The refractive index of a liquid containing dissolved 
gas or gas collected in very small bubbles can be calcu- 
lated using the Bruggman effective medium approxima- 
tion [45[. If the volume fraction of gas is small, fi <C 1, 
then 



An, 



3n (l 



2(1 + 2n 2 ) 



(7) 



where n w 1.34 is the refractive index of the solution at 
A = 633 nm. The volume fraction of gas is expressed as 
fi = Ni/N so i , where iVj is the concentration of z-th gas 
and N so i s» 3.4 x 10 22 cm~ 3 is the concentration of the 
solution. If the gas in the volume cannot disappear, for 
example, due to chemical reaction, then the vibrometer 
signal must be proportional to the total gas flux produced 
by the Faraday current pll ]: 



v(t) = 



dAd 0.35 



dt 



N.. 



[Ji(t) + Mt)] 



(8) 



where v(t) is the signal of the vibrometer and its opposite 
sign is defined by the instrument. The total flux of the 
produced molecules Jtot = Ji + J2 is a non-decreasing 
function of time. 

We observed the process with an instrument Polytec 
MSA-400. The laser beam (A = 633 nm), with a diame- 
ter of 1.5 nm, was focused on the electrode in its center 
at small distance from the edge. The signal is presented 
in Fig. inja) and (b). When frequency of alternating sign 
pulses was low, large bubbles can be formed that scatter 
light significantly. In this case one cannot extract help- 
ful information from the signal as shown in Fig. |9](a) for 
/ = 20 kHz. However, for higher frequencies the bubbles 
small in comparison with X/tt are formed and one can 
clearly see the signal (see FiglHKb)). The signal becomes 
weaker in the inner areas of the electrode in accordance 
with the current distribution. It disappears at all outside 
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FIG. 9: Vibrometer signal, (a) Processing at low frequency, 
/ = 20 kHz, results in formation of large bubbles strongly 
scattering the light. The thick blue curve is the signal and 
the thin red curve is the current, (b) The signal and current 
as functions of time for / = 60 kHz. The peaks in the signal 
are in phase with the voltage (current) pulses, (c) The signal 
in (b) integrated over time (thick blue). The red dashed line 
presents the current and is given here only for eye guidance. 

of the electrodes. The latter demonstrates that the signal 
cannot be related to any kind of vibrations. Integration 
of the signal over time gives the variation of the optical 
path Ad(t). The integrated signal is presented in Fig. 
He). 

The function Ad(t) demonstrates prominent maxi- 
mums that are in phase with the voltage pulses. These 
maximums mean periodic decrease of the gas concentra- 
tion in the liquid. The only reasonable explanation for 
this decrease is the reaction between hydrogen and oxy- 
gen. Of course, the vibrometer signal is sensitive only to 
the overall concentration of gas in the liquid. It is not 
possible to say if the gas interacts as dissolved molecules 
or inside of nanobubbles. However, from observations of 
hydrogen and oxygen on the same time scale we know 
that each gas is collected in nanobubbles. It is hardly 
possible that the gases produced by alternating polarity 
pulses will behave differently. 

The linear trend in Ad(t) means that not all produced 
gas is burned in the reaction but its small part is left 
after each period. This residual gas exists mainly in the 
form of dissolved molecules. When the gas production is 
switched off the residual gas disappears. To all appear- 
ance it disappears also as the result of reaction because 
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FIG. 10: Modification of the electrode surface, (top row) 
The left image shows untreated Pt electrode. The second im- 
age (from left to right) is the same electrode after 55 ms of 
processing. The third and forth images are for W and Au 
electrodes, respectively, after 50 ms of processing, (bottom 
row) Successive images of changes in the Pt electrode for dif- 
ferent processing times. The colormap is changed for better 
visibility. 

reduction of gas concentration due to diffusion has to 
proceed slower than observed. However, this point needs 
more attention. 



C. Modification of the surface of electrodes 

The effect of gas disappearance is accompanied by 
modification of the electrode surface (see Fig. [T0|) . We 
call this phenomenon "wear" of electrodes. It never has 
been observed for single polarity electrolysis but well vis- 
ible for the alternating polarity processing. The effect is 
observed for all investigated materials Pt, Pd, W, and Au 
but manifests itself to a different degree. The strongest 
effect is observed for gold and the weakest one is for tung- 
sten. It is in correlation with the material yield strength. 
The effect increases with the process time as shown in 
the bottom row of Fig. [TTJJ 

The chemical origin of the surface modification can 
be excluded because it is observed for different metals. 
On the other hand, we can exclude the electrochemical 
origin of the modification because it is not observed for 
single polarity pulses. We expect that the surface mod- 
ification happens as the result of local release of energy 
from the exploding nanobubbles that mechanically mod- 
ify and shift the material of electrodes. This assump- 
tion naturally explains correlation with the material yield 
strength. 

We observed only weak modification of the electrodes 
for the alternating pulse electrolysis in the NaCl solution. 
In the KI solution there was no visible modification of 
the electrodes at all. These facts are in good agreement 
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FIG. 11: AFM scan of the modified surface. The process 
time was 55 ms (a) Optical image of the electrodes in air. 
An approximate location of the AFM scan is shown, (b) 3D 
topography of the electrode. Significant displacement of the 
material is observed, (c) Platinum surface before the process 
(normal roughness), (d) Modified roughness of the electrode 
after the process. The roughness changes on both lateral and 
vertical scales. 

with the enthalpy of the reactions between H2 and O2, 
CI2, or I2: the larger the enthalpy the stronger is the 
modification of the surface. 

The assumption of the mechanical modification is sup- 
ported by the analysis of the modified surface with an 
atomic force microscope (AFM). The AFM scan of the 
modified electrode is shown in Fig. 111! One can see 
significant displacement of the material in Fig. [TTT b'l . 
but the most important is that the process changes com- 
pletely the roughness topography of the original platinum 
film as one can see in Fig. ITTT c) and (d). 

The experimental information collected so far on the 
alternating pulse electrolysis (see also [2l|) provides 
strong evidence that the bubbles containing stoichiomet- 
ric mixture of H 2 and O2 gases explode as the result of 
spontaneous reaction between gases if the bubble size is 
smaller than 100 — 200 nm. For pulses of low frequency 
the stoichiometric composition is reached for larger bub- 
bles, which survive and can be observed with the stro- 
boscope. Mechanism of the reaction in small bubbles is 
still not clear. 



V. DISCUSSION 

Let us summarize what we learned about the transient 
nanobubbles using different methods of observation. Of 
course, any bubble growing from a nuclei to its actual size 



gets through the phase when its size is in the nanorange. 
In this paper we have been concentrated on the bubbles 
that grow and live in the nanophase a very short time ~ 
10 /xs. This regime is realized only when the bubbles are 
homogeneously produced nearby the electrode surface. 



A. Evolution of transient nanobubbles 

Homogeneous bubble production is closely related to 
the high current density that results in a huge supersat- 
uration within the diffusion layer. The observed current 
density is large because we consider the system on the 
short-time scale t ~ 10 /its. In this case the current is 
not restricted by the diffusion of electroactive species to 
the surface as it happens for the long-time electrolysis. 

For single polarity pulses the critical nucleus are pro- 
duced by the density fluctuations that is possible be- 
cause of very high supersaturation. These nucleus start 
to grow diffusively. On the time scale r ~ 10 /xs they 
grow to nanobubbles of the size d < 200 nm close to the 
limit of stroboscopic visibility. At longer time the sep- 
arate nanobubbles start to aggregate to form microbub- 
bles. Already at r ~ 100 /xs most of the bubbles exist 
as microbubbles. Small increase of the temperature on 
10 — 20 degrees strongly intensify the bubble nucleation 
and grows as it should be for homogeneous bubble pro- 
duction. At higher temperature the microbubbles ap- 
peared already at t ~ 30 /xs or even earlier. 

When all the gas is collected in microbubbles the su- 
persaturation in the surrounding liquid is low and the 
process becomes very slow. The microbubbles slowly dis- 
solve diffusively. For example, a bubble of 10 /xm in di- 
ameter will dissolve for a time of the order of 100 ms. 

The bubbles produced by the alternating polarity 
pulses evolve differently. The first half of the pulse 
produces high supersaturation of one gas, for example, 
H2 . Nucleus of hydrogen are formed homogeneously and 
grown diffusively. The second half of the pulse produces 
O2 molecules in the same location above the electrode, 
but for oxygen it is energetically preferable to diffuse in 
the existing hydrogen bubble than to form a new O2 
bubble. This is why the alternating polarity pulses pro- 
duce bubbles containing mixture of hydrogen and oxygen. 
This mixture is stoichiometric if positive and negative 
halves of the pulse have the same duration. 

If a hydrogen bubble has grown large enough before 
oxygen starts to diffuse inside, then such a bubble evolves 
similarly to the bubbles containing only one gas. The 
lower the frequency of the alternating pulses the longer 
time has a bubble to grow. The critical frequency is 
around 20 kHz. The critical size of the bubble cannot 
be estimated very precisely from the experiment, but 
roughly it is around d ~ 150 nm 21]. For high fre- 
quency, / > 20 kHz, a bubble has no time to grow large. 
The gas mixture in the bubble becomes stoichiometric 
while its size is still under 150 nm. It was established 
(see also [2l|) that in such a bubble hydrogen and oxy- 
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FIG. 12: (a) Clean electrodes, (b) AFM image of the elec- 
trode surface after fouling with 50 nm nanoparticles. (c) Op- 
tical image of the fouled electrodes, (d) After cleaning with 
homogeneously produced oxygen nanobubbles. Similar effect 
is observed for hydrogen nanobubbles. 

gen react spontaneously with the bubble disintegration 
due to large released energy. 

The mechanism of the reaction inside of nanobubbles is 
not clear. Bulk mixture of gases does not react at room 
temperature and the pressure that is equivalent to the 
Laplace pressure P c = Aj/dtz 30 bar for d — 100 nm. Of 
course, we cannot apply the classical combustion theory 
,46j to a strongly confined system such as nanobubble. 
To all appearance fast microsecond dynamics plays an 
important role for the combustion of mixture of gases. 
The detailed mechanism of the reaction is an important 
point that needs further attention. 

The transient nanobubbles considered in this paper 
cannot help with the understanding of stability of the 
surface or bulk nanobubbles. The stabilization mech- 
anism can work in these bubbles too but they cease to 
exists due to external reasons such as aggregation or com- 
bustion. The transient nanobubbles exist in a very thin 
diffusive layer above the surface, but we cannot say for 
sure if they attached to the surface or not. Relation be- 
tween the transient and surface/bulk nanobubbles is an 
additional point that has to be understood. 

B. Possible applications 

We expect that the transient nanobubbles have a wide 
application potential because of their unique properties. 
For example, dense coverage with homogeneously pro- 
duced nanobubbles can be used for surface cleaning in 
electronic industry and other applications. Figure [T21 
demonstrates the proof of concept, where the original 
clean electrodes were fouled with nanoparticles. The fol- 
lowing homogeneous production of nanobubbles cleans 
the surface. The cleaning happens as the result of aggre- 
gation of nanobubbles that pulls the nanoparticles out 
of the surface. The cleaning effect was already stressed 
for the surface nanobubbles produced electrochemically 
[47| or with other methods [3ll l48l |49| . However, ho- 
mogeneous bubble production is better designed for local 
cleaning and has pronounced dragging effect on nanopar- 
ticles. 

Especially interesting applications are expected for the 
exploding nanobubble phenomenon. Combustion inside 



of a nanobubble results in a significant energy deposition 
that is localized in a nanovolume. It can be applied, for 
example, to delete strongly adhered nanoparticles from 
masks or wafers in the electronic industry. One can easy 
change the bubble size by the electrical means (pulses 
frequency) varying significantly the energy of explosion. 

Combustion in microscopic volumes is considered as 
difficult or impossible due to fast heat escape via the vol- 
ume boundary [501 . The situation becomes different for 
nanovolumes. Exploding nanobubbles are the smallest 
combustion chamber realized in nature [5l| . It opens up 
a principal possibility to fabricate an internal combustion 
engine microscopic in all three dimensions. 



VI. CONCLUSIONS 

In this paper we described a special type of nanobub- 
bles, which we call transient nanobubbles. These bubbles 
exist very short time < 100 fis before aggregating with 
other nanobubbles to form microbubbles or before ex- 
ploding disintegration due to reaction between hydrogen 
and oxygen. The most distinctive feature for all tran- 
sient nanobubbles is their homogeneous production in 
the electrochemical process. In contrast with the nor- 
mal heterogeneous bubbling these bubbles nucleate due 
to fluctuations of liquid density. 

Homogeneous nucleation is possible only for extremely 
high supersaturations that is related to very high current 
densities. These conditions can be met on a shot-time 
scale ~ 10 fj,s when the diffusion of charged particles 
in the electrolyte does not restrict the current. On the 
short-time scale the gas produced on the electrode surface 
has no time to diffuse far away from this surface. For 
this reason high supersaturation can be realized only in 
a diffusion layer of thickness ~ 100 nm. In this thin layer 
only nanobubbles can be formed. 

Short voltage pulses of fixed polarity repeated with 
some frequency / can densely cover the electrode sur- 
face with nanobubbles containing hydrogen or oxygen. 
These nanobubbles aggregate and form microbubbles on 
the time scale ~ 100 (is. The following dynamics of mi- 
crobubbles becomes very slow. 

Short voltage pulses of alternating polarity produce 
nanobubbles containing mixture of hydrogen and oxygen. 
If these bubbles are smaller than 150 nm the gases inside 
of them react spontaneously. As a result of this reaction 
significant energy is released in close proximity of the 
electrode surface. Precise mechanism of the spontaneous 
reaction inside of nanobubbles is still not clear and needs 
clarification. 
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